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Magnetic resonance imaging measures have been proposed as objective markers to study upper motor 
neuron loss in motor neuron disorders. Cross-sectional studies have identified imaging differences between 
groups of healthy controls and patients with amyotrophic lateral sclerosis (ALS) or primary lateral sclerosis 
(PLS) that correlate with disease severity, but it is not known whether imaging measures change as disease 
progresses. Additionally, whether imaging measures change in a similar fashion with disease progression in 
PLS and ALS is unclear. To address these questions, clinical and imaging evaluations were first carried out in a 
prospective cross-sectional study of 23 ALS and 22 PLS patients with similar motor impairment and 19 
age-matched healthy controls. Clinical evaluations consisted of a neurological examination, the ALS Function- 
al rating scale-revised, and measures of finger tapping, gait, and timed speech. Age and ALSFRS score were 
not different, but PLS patients had longer duration of symptoms. Imaging measures examined were cortical 
thickness, regional brain volumes, and diffusion tensor imaging of the corticospinal tract and callosum. Imag- 
ing measures that differed from controls in a cross-sectional vertex-wise analysis were used as regions of in- 
terest for longitudinal analysis, which was carried out in 9 of the ALS patients (interval 1.26 ±0.72 years) and 
12 PLS patients (interval 2.08 ±0.93 years). In the cross-sectional study both groups had areas of cortical 
thinning, which was more extensive in motor regions in PLS patients. At follow-up, clinical measures de- 
clined more in ALS than PLS patients. Cortical thinning and grey matter volume loss of the precentral gyri 
progressed over the follow-up interval. Fractional anisotropy of the corticospinal tracts remained stable, 
but the cross-sectional area declined in ALS patients. Changes in clinical measures correlated with changes 
in precentral cortical thickness and grey matter volume. The rate of cortical thinning was greater in ALS pa- 
tients with shorter disease durations, suggesting that thickness decreases in a non-linear fashion. Thus, cor- 
tical thickness changes are a potential imaging marker for disease progression in individual patients, but the 
magnitude of change likely depends on disease duration and progression rate. Differences between PLS and 
ALS patients in the magnitude of thinning in cross-sectional studies are likely to reflect longer disease dura- 
tion. We conclude that there is an evolution of structural imaging changes with disease progression in motor 
neuron disorders. Some changes, such as diffusion properties of the corticospinal tract, occur early while cor- 
tical thinning and volume loss occur later. 

© 2012 The Authors. Published by Elsevier Inc. All rights reserved. 
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1. Introduction 

Amyotrophic lateral sclerosis (ALS) and primary lateral sclerosis 
(PLS) are clinically distinct motor neuron disorders. Both disorders, 
however, are characterized by the progressive dysfunction and loss of 
corticospinal, or upper, motor neurons (UMNs). Although clinical signs 
are used to diagnose UMN dysfunction, they provide limited information 
about the severity of UMN loss. Several types of tests have been proposed 
as objective measures of UMN dysfunction (Bowser et al., 2006; 
Kaufmann et al., 2004; Turner and Modo, 2010; Ziemann et al., 1997). 
Among these, magnetic resonance imaging (MR1) of the brain has re- 
ceived the most widespread acceptance as a potential method to quanti- 
fy UMN dysfunction (Turner et al., 2011, 2012). Quantitative imaging 
methods, including volumetric measures (Tartaglia et al., 2009; Thivard 
et al., 2007), diffusion tensor imaging (DTI) (Agosta et al., 2010a; Ellis 
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et al., 1999; Iwata et al., 2008, 2011), magnetic resonance spectroscopy 
(Mitsumoto et al., 2007; Pyra et al., 2010), and resting state functional 
MR1 (Agosta et al., 201 1, 2012a; Verstraete et al., 2010), show differences 
between healthy controls and patients with ALS or PLS. In cross-sectional 
studies of ALS patients, some imaging measures correlated with mea- 
sures of disease severity. For example, within a group of ALS patients, 
the fractional anisotropy (FA) of the corticospinal tract in DTI studies 
was correlated with scores on the ALS functional rating scale (ALSFRS-R) 
and with central motor conduction time (Ellis et al., 1999; Iwata et 
al., 2008). From such data it has been suggested that the FA of the 
corticospinal tract could serve as a marker of clinical progression. 

However, it has not been established that quantitative MRI measures 
change with clinical progression in individual patients. Only a small 
number of patients have been studied longitudinally; the relationship be- 
tween imaging findings and clinical measures has largely been inferred 
from cross-sectional studies. The relationship between FA and clinical 
function may be complex. We found that PLS patients had reduced 
corticospinal tract FA, but unlike ALS patients in the same study, measures 
of corticospinal FA were not correlated with measures of UMN function; 
in PLS patients clinical measures correlated better with the mean diffusiv- 
ity (MD) of the corticospinal tract (Iwata et al., 2011). We wondered if 
such differences between PLS and ALS patients were related to differences 
in disease duration or reflected differences in their underlying biology. In 
ALS, the median survival is 3-5 years from the time of diagnosis 
(Czaplinski et al., 2006), whereas the median survival in PLS is often sev- 
eral decades long (Floeter and Mills, 2009; Tartaglia et al., 2007). 

On the other hand, imaging measures could differ between ALS 
and PLS because of differing pathology. A leading hypothesis is that 
both disorders lie on the same disease continuum with a common pa- 
thology characterized by neuronal loss and intra-neuronal inclusions, 
and that the different clinical phenotypes are caused by regional dif- 
ferences in the burden of pathology. Autopsy studies of ALS brains 
have found that neuronal inclusions containing TDP-43 extend to 
extra-motor cortical regions (Geser et al., 2008). This raises the ques- 
tion whether PLS patients have a more restricted phenotype because 
of a more restricted distribution of pathology. Because there are few 
post-mortem studies of PLS patients, it is difficult to test this hypoth- 
esis directly. MRI offers a non-invasive alternative to look for affected 
brain regions of living PLS and ALS patients. To date, only a few imag- 
ing studies have directly compared cohorts of patients with PLS to pa- 
tients with ALS (Ciccarelli et al., 2009; Iwata et al., 2011). 

The first part of the study addressed the hypothesis that affected re- 
gions of the brain differ between PLS and ALS patients, as would be con- 
sistent with a differential distribution of pathology. Imaging measures 
that potentially could reflect pathological changes, such as cortical thick- 
ness and volumes of grey and white matter tracts, were assessed in 
a cross-sectional study comparing brain MRIs from patients to age- 
matched healthy controls. To first visualize areas of cortical thinning, an 
unbiased whole-brain vertex-wise comparison was carried out. Measures 
of the cortical thickness and grey and white volumes of the regions of in- 
terest identified in that unbiased analysis were then compared between 
the patient groups. The second question examined was whether changes 
in imaging measures correlated with clinical progression. Approximately 
half of the patients in the cross-sectional study were able to participate in 
longitudinal imaging studies. In this subset of patients, we looked at the 
correlations between changes in clinical and imaging measures, focusing 
on regions of interest identified in the cross-sectional study and white 
matter tracts previously identified with clinical severity (Agosta et al., 
2010a, 2012b; Ellis et al., 1999; Iwata et al., 201 1 ). 

2. Methods 

2.1. Subjects 

Twenty-two patients who fulfilled Pringle's criteria for PLS (Pringle 
et al., 1992), 23 patients with clinical probable or definite ALS by revised 



El Escorial criteria (Brooks et al., 2000), and 19 age-matched healthy 
controls participated in the cross-sectional study. All subjects gave writ- 
ten, informed consent for the protocol, which was approved by the 
Institutional Review Board. Some of these patients also participated in 
a previously reported diffusion tensor imaging study (Iwata et al., 
201 1 ). A subset of 12 PLS and 9 ALS patients who did not develop con- 
traindications to MRI or significant respiratory or bulbar dysfunction 
had longitudinal clinical evaluations and neuroimaging studies. The in- 
terval between scans was longer for PLS patients than for ALS patients in 
order to allow measureable clinical progression between scans. 

2.2. Clinical evaluation 

A detailed history was obtained to determine the site of onset, 
progression of symptoms, and duration of disease. Patients and care- 
givers were questioned about behavior that would indicate cognitive 
changes or frontotemporal dysfunction. Most patients had neuropsy- 
chological testing. None of the patients had frontotemporal dementia 
or a family history of ALS or frontotemporal dementia. All ALS and PLS 
patients underwent diagnostic testing to confirm that they fulfilled 
the El Escorial criteria for probable or definite ALS (Brooks et al., 
2000) or the Pringle criteria for PLS (Pringle et al., 1992). At each 
visit, a neurologist performed a neurological examination on all sub- 
jects. All healthy controls had normal neurological examinations. 
The revised ALS Functional Rating Scale (ALSFRS-R) was recorded 
for all subjects (Cedarbaum et al., 1999). Finger tapping speed, time 
to walk 20 feet, and time to read a standard passage were measured 
to quantify motor function. At the first evaluation 4 patients were un- 
able to speak, 4 were unable to tap fingers, and 7 patients were un- 
able to walk. For statistical analysis, in patients who were unable or 
became unable to perform these tasks, a value of zero was imputed 
for finger tapping and a value one second greater than that of the 
most severely affected patient was imputed for speech and gait. 
Although the true measures of motor function are unknown, this is 
a conservative method for imputing values. The rate of disease pro- 
gression from onset of symptoms to the first study evaluation was 
estimated by the following calculation: (48 - ALSFRS-R score at first 
clinical evaluation)/time in months between first symptom and first 
study evaluation. The rate of disease progression during the follow- 
up interval was calculated as follows: (ALSFRS-R score at first time 
point evaluation - ALSFRS-R score at last time point)/follow-up inter- 
val in months. 

2.3. Magnetic resonance imaging analysis 

Magnetic resonance imaging (MRI) studies were performed on a 3 T 
scanner (Philips Achieva, Best, the Netherlands) using a receive-only, 
eight-channel SENSE head coil. For volumetric and thickness measure- 
ments, a high resolution Tl -weighted image was obtained using a 
three dimensional turbo field echo sequence (3D Tl TFE) with TR= 
8.6 ms, TE = 3.9 ms, Tl = 700 ms, flip angle = 6°, FOV = 240 mm, matrix 
size = 256x256, slice thickness = 1 mm and 140 axial slices. For DTI 
studies, multi-slice diffusion weighted imaging was acquired with diffu- 
sion weighting along 32 non-collinear directions (b = 1000 s/mm2) and 
one volume without diffusion gradients applied (bO). The sequence was 
repeated 4 times to increase the signal to noise ratio (TE = 86 ms, FOV = 
240 mm, matrix size = 96x96 reconstructed to 128x 128, voxel size = 
1.875 mm x 1.875 mm x 2.5 mm, slice thickness = 2.5 mm, and 55 con- 
tiguous axial slices aligned parallel to the AC-PC line). 

Cortical reconstruction and volumetric segmentation was performed 
with the FreeSurfer image analysis suite (http://surfer.nmr.mgh.harvard. 
edu/). The technical details of these procedures have been described and 
validated in prior publications (Dale et al., 1999; Fischl and Dale, 2000; 
Fischl et al., 1999a, 1999b). Briefly, the processing included skull strip- 
ping, Talairach transformation, optimization of the grey matter-white 
matter and grey matter-CSF boundaries, segmentation, and tessellation 
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(Dale et al., 1999; Fischl and Dale, 2000). The tessellated surfaces were 
then inflated and registered to a spherical atlas which allowed for 
parcellation of the cerebral cortex into units based on gyral and sulcal 
structure, and creation of surface based data, including maps of curvature 
and sulcal depth (Fischl et al., 1999a, 1999b). Cortical thickness was cal- 
culated as the closest distance from the grey/white boundary to the grey/ 
CSF boundary at each vertex on the tessellated surface (Fischl and Dale, 
2000). Cortical gyri were parcellated as described by Desikan et al. 
(2006) for region-of-interest analyses. To measure changes in cortical 
thickness and volume between the initial and follow-up scans, images 
were processed using the longitudinal stream in FreeSurfer (http:// 
surfer.nmr.mgh.hai-vard.edu/fswiki/LongitudinalProcessing) ( Reuter et 
al., 2012). Briefly, this processing steam creates an unbiased within- 
subject template space and image using a robust inverse consistent reg- 
istration (Reuter et al., 2010). Skull stripping, Talairach transformation, 
spherical maps, and parcellation are then initialized from the within- 
subject template for each subject's time point, utilizing common infor- 
mation and thereby increasing reliability and statistical power (Reuter 
et al., 2012). 

Diffusion tensor imaging, and measurements of diffusion proper- 
ties by fiber tracking, were carried out using methods previously 
described (Danielian et al., 2010; Iwata et al., 2011). Initial image 
processing, including eddy current correction, linear registration, 
and masking, was performed using FSL (University of Oxford, UK) 
(Smith et al., 2004). The diffusion tensor was calculated and fiber 
tracking was done using DtiStudio (Jiang et al., 2006) (www.MriStudio. 
org) using the Fiber Assignment by Continuous Tracking (FACT) method 
(Mori et al., 1999; Xue et al., 1999) with a minimum FA value of 0.2 for 
tracking termination and maximum 40 degree angle between two adja- 
cent eigenvectors. The corticospinal tract (CST) was tracked using a mul- 
tiple region-of-interest (ROI) approach (Iwata et al., 201 1 ; Wakana et al., 
2007) previously shown to have good reliability for longitudinal studies 
(Danielian et al., 2010) and the corpus callosum (CC) was segmented 
into the motor fibers of the corpus callosum, genu and splenium 
according to the technique of Hofer and Frahm (2006). Whole-tract 
measures of mean fractional anisotropy (FA) and mean diffusivity 
(MD; 10~ 3 mm 2 /s) were calculated within each delineated tract. To as- 
sess for changes in cross-sectional area along the rostro-caudal extent of 
the corticospinal tract, profiles of the mean cross-sectional area within 
the tracked CST in each axial slice were constructed as previously de- 
scribed (Iwata et al., 2011). 

In addition, longitudinal changes in diffusion metrics of the CST in 
the same anatomic space were evaluated by an atlas-based method 
(Oishi et al., 2009). DTI data was processed with DiffeoMap and 
RoiEditor software (Jiang et al., 2006) (www.MriStudio.org). The im- 
ages were normalized to the ICBM-DTI-81 coordinates (Mori et al., 2008) 
using 12-mode affine transformation. Subsequently, a non-linear trans- 
formation using dual-contrast large deformation diffeomorphic metric 
mapping (LDDMM) was applied to register the FA and Trace images 
to the 1 mm isotropic resolution template (181x217x181 matrix). 
Because both linear and non-linear transformations are reciprocal 
procedures, the inversely transformed WMPMI parcellation map 
was superimposed onto the original DTI metrics maps from each 
patient, leading to the parcellation of the brain into 176 anatomical 
structures. For segmentation, a threshold of FA<0.2 and Trace > 0.0045 
was applied to exclude grey matter and CSF from white matter on the 
corresponding images. FA and MD measures were analyzed in five 
anatomically defined regions of the CST parcellated by this atlas-based 
method, shown in Fig. 1. 

2.4. Statistical analysis 

All data in tables are presented as mean ± standard deviation. 
Basic descriptive statistics and statistical analyses were performed 
using SPSS Statistics version 17.0 or SAS version 9.2. For clinical mea- 
sures, Box-Cox transformation was applied to timed speech (natural 



logarithm), gait measurements (inverse square root) and ALSFRS-R 
(square); Shapiro- Wilk was used to test normality (based on the 
residuals). Demographic and clinical data were compared between 
groups in the cross-sectional study by one-way unequal variance 
ANOVA with post hoc Tukey's test. In the longitudinal study, within- 
subject changes in clinical measures were determined using a mixed 
model, described below. 

2.4.1. Cross-sectional imaging analysis 

For the cross-sectional study, unbiased whole-brain statistical 
comparisons were performed with FreeSurfer's QDEC application 
(www.surfer.nmr.mgh.harvard.edu). Surface maps depicting regions 
with statistically significant differences in cortical thickness at each 
vertex across the entire cortical mantle were determined in QDEC 
with general linear models (GLMs) using p<0.001 uncorrected as 
the threshold for a significant cluster. Pain/vise comparisons were 
made between ALS patients, PLS patients and healthy controls. Age 
and gender were not included as covariates in the QDEC model be- 
cause these measures did not differ between the groups. Volumes of 
grey and white matter regions in which thinning was identified by 
QDEC were compared with a one-way ANOVA with post hoc Tukey's 
test. 

2.5. Longitudinal analysis 

Cortical regions found to have reduced thickness in the cross- 
sectional QDEC analysis were used as regions-of-interest in the longi- 
tudinal subset. Cortical thickness, grey matter volume and white mat- 
ter volume of those regions were assessed, as well as changes in the 
FA and MD of the corticospinal tract and mid-posterior portion of 
the corpus callosum. The change between the two visits in these im- 
aging measures and clinical measures was tested by a mixed model 
with patient as random effect. The full model included visits, diagno- 
sis, interaction of diagnosis by visit, age and scan interval (p-value of 
0.1 used as model selection). For those clinical and imaging measures 
showing significant change between the two visits, a general linear 
model was applied to examine the within-subject correlation, which 
described an increase in a clinical measure within the individual 
that was associated with an increase (or decrease) in an imaging 
measure by removing the difference between subjects (Bland and 
Airman, 1995). When an interaction between a clinical measure and 
diagnosis was significant, the association analysis was applied to 
ALS and PLS groups separately. 

Brains processed with the FreeSurfer longitudinal stream were 
analyzed in QDEC with diagnosis as a discrete factor and disease 
duration as the covariate to generate surface maps showing where 
rates of thinning (thickness2 — thicknessl ) / (time2 — timel ) differed 
between diagnoses. 

3. Results 

3.1. Cross-sectional demographics and clinical testing 

There were no differences in age or gender between the three 
groups in the cross-sectional study (Table 1, left panel). The ALSFRS-R, 
timed speech, gait, and finger tapping rate were significantly different 
between healthy controls and patients (p<0.001), but did not differ 
between the two patient groups. As expected, the mean duration of 
disease was longer in the PLS group than in the ALS group and the 
rate of disease progression was less in PLS patients than ALS patients. 
In the subset of 21 patients who underwent longitudinal studies, the 
age and gender ratios were similar to the cross-sectional group, but 
there was a trend towards slightly better clinical measures and a slower 
initial rate of disease progression in this subset compared to the 
cross-sectional group as a whole (Table 1, right panel). 
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Fig. 1. Regions of the corticospinal tract defined by an atlas-based method that were used for longitudinal analysis of diffusion tensor imaging. Labeled regions include the pons 
(CST, blue), cerebral peduncle (Peduncle, turquoise), posterior limb of the internal capsule (PL1C, green), superior corona radiata (SCR, yellow), and the subcortical pre-central 
white matter (PreCG, red) (Oishi et al., 2009). 



3.2. Cross-sectional imaging 

3.2.1. Group differences in cortical thickness and grey and white matter 
volumes 

Maps depicting pair-wise regional differences in cortical thickness 
between patients and healthy controls are shown in Fig. 2. In PLS pa- 
tients cortical thinning was mostly limited to motor regions, with five 
clusters of thinning in the left precentral cortex and two clusters of 
thinning in the right precentral cortex (Fig. 2B). The coordinates of 
the maxima of the clusters in the precentral region corresponded ap- 
proximately to Brodmann area 4 in Talairach space (Table 2) (Fischl 
et al., 2008). There was only one very small cluster of thinning outside 
of the motor area in PLS, in the left lateral orbital frontal region. In the 
ALS patients, clusters of cortical thinning were less extensive than in 
PLS patients (Fig. 2A). There were two clusters of cortical thinning in 
the left and right motor areas, with coordinates corresponding to 
Brodmann areas 4 and S in Talairach space (Table 2). Two clusters 
of cortical thinning occurred outside motor areas in ALS, in the left 
postcentral cortex and a small cluster in the rostral middle frontal 
cortex (Table 2). Based on these locations of thinning in the whole- 
brain maps, the following seven parcellated regions were selected 
for further region-of-interest analyses: right and left precentral 
gyrus, right and left paracentral region, left postcentral gyrus, left ros- 
tral middle frontal region, and left lateral orbital region. 

The total brain, white, and grey matter volumes did not differ be- 
tween groups ( ANOVA, p > 0.2 ), indicating that there was no global at- 
rophy. However, the grey and white matter volumes of some regions of 



interest, normalized for intracranial volume, differed between groups, 
suggesting focal atrophy (Supplementary Table 1). The right and left 
precentral grey matter volume was reduced in PLS and ALS patients 
and the grey matter of the post-central gyrus was reduced in ALS pa- 
tients compared to controls. PLS patients also had reduced right and 
left precentral and paracentral white matter volume and ALS patients 
had reduced volume of the postcentral gyrus white matter. 



3.3. Longitudinal study 

3.3.1. Longitudinal changes in clinical and motor measures 

Longitudinal evaluations were carried out in the subset of 9 ALS 
patients and 12 PLS patients who did not develop contraindications 
to MRI or significant respiratory dysfunction. The mean interval 
between follow-up evaluations was 1.26 ±0.72 years for the ALS 
patients, and 2.08 ± 0.93 years for the PLS patients. Both groups exhibited 
a decline in the ALSFRS-R scores and in several of the motor measures 
over the follow-up interval, although the decline was more severe in 
the ALS patient group (Table 1, right panel). Mean ALSFRS-R scores de- 
clined by approximately 6 points in the ALS patient group and by 3.5 
points in the PLS patient group. However, this includes five PLS patients 
who had stable ALSFRS-R scores that dropped one point or less. Measures 
of finger tapping speed declined bilaterally significantly in ALS patients 
over the follow-up interval. Two ALS patients became unable to perform 
one or more of the motor measures at follow-up that were initially 
performed. In PLS patients, only left-hand finger tapping exhibited 
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Table 1 

Characteristics of subjects in cross-sectional study and longitudinal subset. 





Cross-sectional study 




Longitudinal subset 
PLS 


ALS 






Control 


PLS 


ALS 


Initial 


Follow-up 


Initial 


Pf\\ 1 n\Aj— i i n 
ruiiuvv uu 


Male:Female 


11:8 


8:14 


12:11 


7:5 




5:4 




Age (years) 


58.7 ±5.9 


59.0 ±7.9 


55.8 ±11.3 


56.9 ±9.4 


59.41 ±9.6 


57.2 ±12.6 


58.38 ±13.1 


Disease duration (years) 




a 13.5±6.7 


3.0 ±2.5 


9.39 ±4.2 


11.87 ±4.3 


3.66 ±2.9 


4.87 ±3.1 


ALSFRS-R 


b 48.0±0.0 


36.3 ±5.2 


33.7±7.1 


39.8 ±3.3 


c 36.2±5.3 


40.2 ±6.3 


c 34.1±9.8 


Speech (s) 


b 51.1±9.0 


79.3 ±26.9 


84.3 ±32.3 


78.8 ±26.0 


81.8 ±26.7 


69.4 ±3 1.7 


85.2 ±41 .8 


Gait (s) 


b 4.2±0.7 


21.4 ±24.0 


44.9 ±49.5 


12.4 ±8.0 


14.8±11.1 


24.9 ±40.8 


54.9 ±57.1 


Finger tapping 


b 86.7±12.7 


48.4 ±14.0 


45.0 ±32.5 


51.0±15.5 


48.4 ±13.3 


63.0±31.8 


c 23.0±28.9 


Rhand (taps/15 s) 
















Finger tapping 


b 85.1±11.6 


41.8 ±14.2 


43.4 ±28.1 


46.4 ±15.0 


c 41.7±13.00 


57.4 ±29.3 


c 23.2±29.3 


Lhand (taps/15 s) 
Progression rate* 


















a 0.09±0.06 


0.77 ±1.0 


0.08 ±0.04 


0.14±0.19 


0.59 ±1.26 


0.51 ±0.58 



Follow-up progression rate = (ALSFRS-R1 — ALSFRS-R2) /follow-up interval in months. 
a PLS and ALS significantly different, p<0.05. 

b Controls significantly different from patients, p<0.001, no difference between PLS and ALS. 

c Significant longitudinal change, p<0.05. 

* Progression rate = (48-ALSFRS-R)/Disease duration in months. 



significant, though modest, declines. There was a trend toward slower the FA and MD remained stable during the follow-up interval in both 
gait in both groups (ALS p = 0.053; PLS p = 0.047). patient groups (not shown). 



3.3.2. Longitudinal changes in diffusion tensor imaging 

The corticospinal tract was analyzed both by a fiber tracking 
method and an atlas-based method. Fiber tracking showed a reduced 
mean FA and increased MD of the corticospinal tract compared to 
previously reported controls (Iwata et al., 2011), but the mean CST 
FA and MD remained relatively stable over the follow-up period, as 
illustrated for individual ALS and PLS patients (Supplementary Fig. 
1 ). However, the mean cross-sectional area of the tracked CST became 
thinner along its rostro-caudal extent in the ALS patient group be- 
tween the initial and follow-up evaluations (Fig. 3A). In PLS patients, 
the cross-sectional area of the CST remained stable over time (Fig. 3B). 
The atlas-based method likewise found no consistent pattern of change 
in the FA and MD measurements over the follow-up interval for the five 
regions of the corticospinal tract defined anatomically (Supplementary 
Fig. 2). The mid-portion of the corpus callosum, where fibers connecting 
the motor cortices travel (Hofer and Frahm, 2006) had lower FA and 
higher MD than previously published controls (Iwata et al., 2011 ), but 



3.3.3. Longitudinal changes in cortical thickness 

Imaging measures from the two time points were obtained using the 
longitudinal processing stream in FreeSurfer (Reuter et al., 2012) and 
analyzed in SAS 9.2. Changes in the mean cortical thickness for each of 
the seven regions of interest (Table 2) were evaluated using a mixed 
model that included factors of age, scan interval, and diagnosis. PLS 
and ALS patients both had significant thinning of the right and left 
precentral gyrus between the two scans (p<0.001 ). There was a signif- 
icant effect of the scan interval (p< 0.005), with a greater degree of thin- 
ning when the interval between scans was longer. Significant thinning 
between scans also occurred in the left postcentral gyrus (p = 0.0095) 
where thinning was affected by age (p = 0.0194), but not by scan inter- 
val. Age did not significantly affect the thinning of the other regions 
studied. 

To assess the effects of disease duration on the rate of thinning in 
ALS and PLS patients, a whole-brain vertex-wise analysis was carried 
out in QDEC using p< 0.001 as the threshold for significance. The 




Fig. 2. Areas of cortical thinning identified in the cross-sectional study by a whole-brain vertex-wise analysis, shown projected on a smoothed average of all brains in this study. 
(FreeSurfer QDEC) Regions of thinning compared to healthy controls are shown for (A) ALS patients and (B) PLS patients. The color of the scale bar indicates the statistical signif- 
icance of the cluster, ranging from the threshold value of p<0.001 to p<0.00001. White lines indicate the approximate region parcellated as precentral gyrus and the yellow lines 
indicate the region parcellated as paracentral region. Hemispheres are viewed from left to right as: left lateral, left medial, right medial, right lateral. 
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Table 2 

Size and coordinates of clusters of cortical thinning in patients compared to controls 
identified in cross-sectional study vertex-wise analysis. 





Cluster size 
(mm 2 ) 


r 

Coordinate ot maxima 
(Talairach coordinates) 


rreeSurter Region 


PLS<Controls 






Left 


273.4 a 


-32.8 -16.3 43.9 


Precentral 




72.6 a 


-50.7 -5.4 41.7 


Precentral 




23.4 a 


-57.6 2.2 22.3 


Precentral 




15.3 a 


-25.9 -11.1 50.1 


Precentral 




8.1 a 


-28.7 -21.4 49.3 


Precentral 




1.8 a 


-30.5 28.4 -15.2 


Lateral orbital frontal 


Right 


83.2 a 


J.D — 1 D.o Do. 1 


Precentral 




54.4 a 


36.3 -10.8 52.9 


Precentral 


ALS<Controls 






Left 


64.5 


-37.9 -10.3 46.8 


Precentral 




40.7 


-7.5 -21.6 69.0 


Paracentral 




1.4 


-39.9 49.6 -8.8 


Rostral middle frontal 




84.0 


-24.5 -26.9 53.8 


Postcentral 


Right 


38.6 


38.6 -9.2 55.5 


Precentral 




6.1 


5.4 -18.4 69.8 


Paracentral 


ALS>Controls 






Right 


43.5 


18.614.2 -19.3 


Lateral orbital frontal 



Threshold for significance was p<0.001 (uncorrected). 
a Comparison also significant at p<0.05 with FDR correction for multiple comparisons. 



resulting cortical maps showed five small areas where the average 
rates of thinning differed between ALS and PLS patients (Fig. 4A). 
These areas were located in the left and right precentral gyrus, the 
right paracentral region, the right superior frontal region, and the 
left cuneus region. For each of the five areas, the slope of the rate of 
thinning versus time from disease onset (using the mid-interval 
time point) differed for ALS patients and PLS patients. In ALS patients 
the rates of thinning in these areas was greatest in patients with a 




Pons 



Cortex 



Fig. 3. Longitudinal changes in the cross-sectional area of the corticospinal tract Dashed 
lines indicate the initial scan visit and solid colored lines indicate the follow-up. Fiber track- 
ing was used to define the CST from the level of the mid-pons to the motor cortex and the 
two sides were averaged for the profiles. (A) In ALS patients, the cross-sectional area of the 
CST declined between the initial and follow-up evaluation. (B) In PLS patients, the CST 
cross-sectional area is smaller than in ALS, but there is no change during the follow-up in- 
terval. Dashed lines from left to right indicate three anatomic levels of the corticospinal 
tract: peduncle, internal capsule, and subcortical white matter. 



shorter disease duration, whereas the rates of thinning in PLS patients 
were similarly low in patients across the range of disease durations 
(Fig. 4B). Thickness of the right and left precentral cortex at each 
time point is shown for individual patients in Fig. 5 to further illus- 
trate these differences. 

3.3.4. Longitudinal changes in volumes of grey and white matter in 
regions of interest 

In most of the regions with significant changes in cortical thick- 
ness, there generally was reduced volume of grey and white matter 
(Supplementary Table 2). The volume of grey matter declined in 
the precentral gyri bilaterally and the left postcentral gyrus in both 
ALS and PLS patients. Both groups also had loss of white matter 
volume in the left postcentral gyrus. PLS patients had loss of white 
matter volume in the precentral gyrus. There was also volume loss 
in the midposterior portion of the corpus callosum, adjusted for 
age. Age was not a significant factor in the volume loss observed in 
the other grey and white matter regions. 

3.3.5. Correlation between longitudinal changes in clinical and imaging 
measures 

Correlations were calculated between clinical and imaging mea- 
sures that had significant declines over the follow-up interval. The 
decline in the ALSFRS-R score was correlated with thinning of the 
precentral cortex (left r= 0.459, p = 0.036; right r = 0.542, p = 0.011), 
loss of precentral grey matter volume (left r = 0.493, p = 0.023; 
right r = 0.560, p = 0.011), loss of left precentral white matter vol- 
ume (r = 0.569, p = .007), and loss of left postcentral white matter 
volume (r = 0.628, p = 0.002) in ALS and PLS patients. In PLS pa- 
tients, the decline in left hand finger tapping speed was associated 
with a decline in right precentral cortex thickness (r = 0.586, p = 0.036) 
and grey matter volume (r= 0.644, p = 0.0176), with a trend toward an 
association with volume of the mid-posterior region of the corpus 
callosum (r= 0.522, p = 0.067). In ALS patients, the decline in right 
hand finger tapping rate was associated with a decline in the white mat- 
ter of the postcentral gyms (r = 0.713, p = 0.0207). 

4. Discussion 

There is a longstanding debate whether primary lateral sclerosis and 
amyotrophic lateral sclerosis are variants of the same disease (Rowland, 
1999; Swash et al., 1999). Two characteristics that distinguish PLS from 
ALS are long survival and lack of significant lower motor neuron signs. 
Both disorders, however, have progressive upper motor neuron loss. 
In ALS, neuronal loss and pathological inclusions are widespread in 
the brain, but it is not known whether the same is true in PLS. The 
first part of this study compared the distribution of cortical thinning, a 
potential imaging marker of pathological change, in ALS and PLS pa- 
tients to healthy controls in a whole-brain analysis. In the motor re- 
gions, cortical thinning occurred in a similar distribution in ALS and 
PLS patients, but the magnitude of thinning was greater in PLS patients. 
In PLS patients, cortical thinning was highly localized to the precentral 
gyrus and adjacent paracentral region, with only one very small cluster 
in the lateral orbitofrontal cortex. Thinning of the motor cortex is con- 
sistent with our earlier report in PLS patients using manual measure- 
ments of the precentral gyms from Tl weighted MRI scans (Butman 
and Floeter, 2007), but that study did not examine the whole brain to 
look for changes in other regions. The current study did not find global 
atrophy in PLS patients, in contrast to a prior report (Tartaglia et al., 
2009). However, in regions with large areas of thinning, i.e. the 
precentral cortex, the volumes of grey and white matter were less 
than in controls, indicating focal atrophy. Our findings do not sup- 
port the hypothesis that cortical pathology is widespread in PLS pa- 
tients. Despite their long duration of disease, imaging differences 
were mostly restricted to the motor regions in PLS and did not spread 
to adjacent cortical regions. Small areas of cortical thinning within 
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Fig. 4. Differences in rates of cortical thinning as a function of disease duration in PLS and ALS patients in a whole-brain vertex-wise analysis. (A) Regions where the average rate of 
change in thickness differs between ALS patients and PLS patients are shown on the smoothed average brain template (FreeSurfer QDEC). The color of the scale bar indicates the 
statistical significance of the cluster, ranging from the threshold value ofp<0.001 top<0.00001. Hemispheres are viewed from left to right as: left lateral, left medial, right medial, 
right lateral. (B) Plots show the relationship between the rate of change in thickness (mm/year) for ALS patients (red squares) and PLS patients (blue circles) and the time of symp- 
tom onset (disease duration at the mid-point of the interval between scans). Data from each plot corresponds to the region of significance that is circled in the brain image directly 
above in (A). ALS patients (red dashed lines) had a faster average rate of thinning for the circled region compared to PLS patients (blue dashed lines). The solid lines show the fitted 
slopes for each disease group. Rate of change in thickness (mm/years) = (thickness2 — thicknessl)/(time2 — timel). 



the precentral gyrus were seen in the ALS patient group, in agreement 
with other studies (Agosta et al„ 2012b; Roccatagliata et al., 2009; 
Verstraete et al., 2012). Thinning was also seen in the postcentral 
gyrus and a small region in the rostral frontal cortex in ALS. In the 
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Fig. 5. Longitudinal changes in the mean cortical thickness of the precentral gyrus in 
patients with ALS (red squares) and PLS (blue circles). Time is referenced to the onset 
of the patient's symptoms. (A) Left precentral gyrus. (B) Right precentral gyrus. 



precentral cortex and the postcentral gyrus, where clusters of thinning 
were largest in ALS, the volume was also reduced. Although thinning 
extended beyond the motor cortex in ALS, it was not as widespread as 
has been described in other reports (Agosta et al., 2012b; Verstraete 
et al., 2012). Although some of the differences between the studies 
can be attributed to differences in the statistical models and method 
for selecting regions of interest, another contributing factor may be 
that none of the ALS patients in the current study had dementia, which 
is associated with greater frontotemporal pathology (Strong and Yang, 
2011). 

The longitudinal portion of the study examined whether imaging 
measures change with clinical progression, focusing on the regions 
identified in the cross-sectional study. In the subset of patients who 
had longitudinal evaluations, cortical thinning was found to progress 
over time. After adjusting for covariates, in both ALS and PLS patients, 
only the precentral gyri showed significant interval change in thick- 
ness. The decline in thickness of the precentral cortex was correlated 
with a decline in the ALSFRS-R score. In PLS patients it also correlated 
with a decline in finger tapping speed. Thinning of the precentral cor- 
tex progressed at a faster rate in ALS patients than PLS patients, and 
within the group of ALS patients, more rapid thinning occurred in 
those with a shorter disease duration. In contrast, in PLS patients 
the rate of thinning was low and relatively similar across the range 
of disease durations. The different rates of thinning could reflect differ- 
ences in the underlying disease processes in ALS and PLS. However, an- 
other interpretation of these data would be that thinning occurs in a 
non-linear fashion, with an early rapid phase of decline followed by a 
longer, slower phase of thinning, in a similar manner in both ALS and 
PLS patients. The detection of significant changes in thickness would 
therefore depend on the timing of the scans in relation to the disease 
duration and the patient's progression rate. Of note, however, is that 
patients with slower rates of progression may be over-represented in 
longitudinal studies, because they are more likely to tolerate repeat 
evaluations and are less likely to drop out than patients who have 
rapid progression of disease. This selection bias was true in the current 
study, as evidenced by the slower progression rates in the subset of 
patients who had follow-up imaging compared the whole group. 
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In the diffusion tensor imaging of the corticospinal tract, longitudi- 
nal changes occurred in cross-sectional area rather than in white matter 
diffusion properties. In this study, we analyzed the corticospinal tract by 
tractography and by an automated image parcellation/atlas-based anal- 
ysis. Both methods yielded the same results, that changes in measures 
of corticospinal FA and MD did not correlate with changes in clinical 
measures. Tractography, or fiber tracking, depends on sufficient axonal 
integrity to produce directionality of diffusion for reconstructing tracts. 
Although fiber tracking of the corticospinal tract had high reliability in 
longitudinal studies of healthy controls (Danielian et al., 2010), it has 
the drawback that in patients with axonal disruption and loss, it can 
only assess surviving fibers. The fiber tracking method found that 
whole-tract measures of corticospinal tract diffusion properties were 
stable during the follow-up period in ALS and PLS patients, although 
they differed from control values at both time points. In ALS, the 
cross-sectional areas of the corticospinal tract declined between scans, 
presumably due to axonal loss and dropout of fibers from tracking. 
The atlas-based method is complementary to the fiber tracking method, 
as it assesses the brain in a structure-by-structure (group of pixels) 
basis using a common panellation map. The method uses large defor- 
mation diffeomorphic metric mapping (LDDMM) as the normalization 
algorithm. The LDDMM has a high accuracy in matching the shapes of 
different brains to a standard atlas. The measurements acquired with 
the atlas-based analysis derive directly from the native images of 
each individual, reducing the imprecision produced by the imaging 
postprocessing. The atlas-based method showed no change in FA and 
MD of the corticospinal tract at the follow-up scan. Although the two 
methods did not find an interval change in diffusion properties, the 
patients' CST FA and MD differed from controls, indicating that changes 
in diffusion properties of the white matter tracts had already occurred 
prior to the time of the first scan. This finding that diffusion properties 
remained relatively stable differs from a study that found a decline in 
FA in the superior portion of the right corticospinal tract with a 
tractography-based region-of-interest analysis which divided each 
CST into three portions (Zhang et al., 2011). In that study, the changes 
were driven by a subset of patients with left-onset ALS (Zhang et al., 
2011). This region is not well-tracked with our tractography method, 
which used a slightly higher FA threshold for truncation of tracking. 
Nevertheless, the atlas-based method also did not find changes in 
these regions. It is possible that the methods we used were not sensitive 
to small changes, particularly given the small sample size in the longitu- 
dinal study. Another possibility is that changes in CST diffusion proper- 
ties occur very early in ALS, and thus were captured in the study by 
Zhang et al. (2011) where the mean disease duration of the group 
was 20.8 months, but not in our study where ALS patients had a mean 
disease duration of 36.6 months at the time of the first scan. 

The tissue changes that give rise to imaging changes such as cortical 
thinning result from a remodeling process that requires time. In addi- 
tion to shrinkage and loss of cortical neurons that could contribute to 
cortical thinning (Hudson et al., 1993; Kiernan and Hudson, 1991, 
1993), there are also inflammatory changes with T-cell infiltration 
(Henkel et al., 2009), reactive astrocytosis (Murayama et al., 1991; 
Vargas and Johnson, 2010), microglia invasion (Turner et al., 2004), 
and iron accumulation (Kwan et al., 2012) in the cortex. These may con- 
tribute to a gradual thinning that continues beyond the period of neuro- 
nal loss, by clearing of necrotic debris and by maturation of glial scarring 
in patients who survive longer. This sequence of changes, and the time 
required for them to occur, would be consistent with the greater extent 
and degree of thinning of the motor cortex in PLS patients compared to 
ALS patients. A similar evolution of tissue changes in the white matter 
could also explain findings from diffusion tensor imaging studies. Nearly 
all cross-sectional diffusion tensor imaging studies have reported 
reduced fractional anisotropy of the corticospinal tract in ALS and PLS 
patients compared to healthy controls, and many studies found a corre- 
lation between the CST FA and clinical severity scores or progression 
rates (Ciccarelli et al., 2006, 2009; Ellis et al., 1999; Iwata et al., 2008, 



2011; Mitsumoto et al., 2007; Sage et al., 2007; Ulug et al., 2004; 
Unrath et al., 2010). In the few longitudinal studies in ALS, lower 
corticospinal tract FA correlated with faster rates of disease progres- 
sion or shorter survival, where ALS patients were followed for one 
(Nickerson et al., 2009) or up to approximately 3.5 years (Agosta et 
al., 2010b). The relatively short follow-up periods that can be achieved 
with ALS patients may not be sufficient to observe the full evolution of 
white matter changes that are seen in PLS patients, some of whom 
have decades-long disease durations. The changes in fractional anisot- 
ropy may be a relatively early event in motor neuron degeneration, 
during active degeneration of corticospinal axons, with loss of cross- 
sectional area of the tract occurring afterward. Others have also noted 
that diffusion properties can change independently of volumetric 
changes (Canu et al., 2011). 

Several cross-sectional studies have reported differences in diffu- 
sion properties in motor fibers of the corpus callosum in ALS and 
PLS compared to healthy controls (Ciccarelli et al., 2009; Douaud et 
al., 201 1 ; Iwata et al., 201 1 ). It has been suggested that reduced frac- 
tional anisotropy in the mid-posterior corpus callosum is specific for 
ALS, and may precede the onset of clinical signs of upper motor neu- 
ron dysfunction (Filippini et al., 2010; Sach et al., 2004). This study 
confirms and extends that proposal, finding that the diffusion proper- 
ties of the callosal motor fibers remained stable over the follow-up in- 
terval, although volume loss occurred. Thus the sequence of imaging 
changes in callosal fibers in motor neuron disease is similar to the 
corticospinal tract, with early changes in diffusion properties and 
later changes in volume. 

This study's finding that progressive thinning and volume loss oc- 
curs primarily the precentral cortex is consistent with a longitudinal 
morphometry study that showed atrophy of the precentral region 
over a 9-month span in ALS patients (Agosta et al., 2009), but differs 
from a recent longitudinal study of ALS patients that found progres- 
sive thinning of extra-motor regions over a 3-10 month follow-up 
interval (Verstraete et al., 2012). One of the limitations of our longitu- 
dinal study is that the small sample sizes limited the number of 
regions that could be assessed in our statistical model. Regions of in- 
terest were selected from areas identified in the cross-sectional 
model which differed from age-matched controls, and did not include 
the whole brain. Of the three extra-motor regions assessed in our lon- 
gitudinal model, only the postcentral gyrus showed progressive thin- 
ning in ALS patients. Several studies have reported more widespread 
atrophy in motor neuron disease (Grosskreutz et al., 2006; Tartaglia 
et al., 2009; Turner et al., 2007), including volume loss in the frontal 
and temporal lobes in patients with cognitive impairment (Chang et 
al., 2005; Strong et al., 2003). In our cross-sectional whole-brain anal- 
ysis, a small cluster of thinning was seen in the lateral orbital frontal 
cortex in PLS patients, and small cluster of thinning occurred in the 
frontal cortex in ALS patients. The lack of more extensive findings in 
extramotor areas may reflect better cognitive function, as there is a 
selection bias in our study for patients who were able to travel to par- 
ticipate in research, compared to a more population-based cohort 
such as from a regional clinic (Grace et al., 201 1 ; Murphy et al., 2008). 

5. Conclusion 

The findings from this longitudinal study favor the idea that there is 
a temporal evolution of imaging changes in the brain in motor neuron 
disorders as upper motor neuron degeneration progresses. We propose 
that alterations in the diffusion properties of the corticospinal tract and 
motor region of the corpus callosum are among the early changes, 
which eventually reach a nadir. Cortical thickness declines at a rapid 
rate initially, then at a slower rate, with progressive atrophy of the af- 
fected grey and white matter regions. Within this framework, longitudi- 
nal changes in imaging measures will be dependent on disease duration 
and progression rate and differences between ALS and PLS patients may 
be explained by differences in disease duration. Future studies with 
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larger sample sizes and additional longitudinal evaluations will be 
needed to test this proposed sequence of imaging changes. In the 
cross-sectional portion of this study, we found small differences be- 
tween PLS and ALS patients in the distribution of cortical thinning, not 
a substantial differential distribution of pathology. However, this con- 
clusion is limited by the clinical similarity of the two patient groups, 
with less widespread pathology in the ALS group than reported in 
other studies. A more diverse ALS patient sample should be studied to 
better assess regional differences between PLS and ALS. Lastly, we con- 
clude that the imaging measure that best correlates with clinical pro- 
gression is the decline in cortical thickness of the precentral gyrus. 

Supplementary data to this article can be found online at http:// 
dx.doi.org/10.1016/j.nicl.2012.12.003. 
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